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SECTION  I 


INTRODUCTION 

During  1965,  deep-well  seismic  records  were  obtained 
from  the  Trigg  well  near  Grapevine,  Texas,  (Figure  1)  and  the  Carter 
well  at  the  Uinta  Basin  Seismological  Observatory  near  Vernal,  Utah, 
(Figure  2).  Data  of  interest  to  the  present  study  are  discucsed  in  this 
special  report. 

At  Grapevine,  the  array  consisted  of  a  vertical-component 
seismometer  at  the  surface;  six  vertical-component  seismometers  locked 
into  the  well  at  depths  of  3500,  4500,  5500,  6500,  7500,  and  8500  ft;  and 
two  horizontal -component  seismometers  at  the  surface.  High-gain  re¬ 
cordings  of  the  above  instruments  were  made  using  channels  1  through 
9,  respectively,  and  simultaneous  low-gain  recordings  were  made  with 
channels  10  through  18,  respectively. 

At  UBO,  outputs  were  recorded  from  a  vertical-component 
seismometer  at  the  surface,  a  north-south  seismometer  at  the  surface, 
an  east-west  seismometer  at  the  surface,  and  six  vertical-component 
seismometers  locked  into  the  well  at  depths  of  3900,  4900,  5900,  6900, 
7900,  and  8900  ft.  The  outputs  of  these  nine  instruments  were  recorded 
in  channels  21  through  29,  respectively.  Low-gain  recordings  were 
made  in  order  to  obtain  teleseismic  records  which  did  not  overmodulate. 
High-gain  recordings  were  made  to  produce  noise  samples  with  maximum 
dynamic  range. 
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Both  *  “cording  programs  used  a  modified  Tex.s  Instruments 
Digital  Field  System  (DFS)*  to  record  multichannel  data  on  magnetic  tape  in 
Texas  Instruments  Automatic  Computer  (TIAC)*  format  with  a  sampling 
interval  of  24  msec.  Grapevine  amplifier  settings  were  such  that,  within 
each  set  of  field  data,  outputs  of  the  deep-well  seismometers  were  subjected 
to  a  gain  factor  of  JO  <i.  e.  .  -20  db)  relative  to  the  output  of  the  vertical- 
component  seismometer  at  the  surface.  Data  were  recorded  in  the  field  in 
30-min  records.  Editing  on  the  TIAC  computer  produced  a  set  of  2-min 
(Grapevine)  and  4-min  (UBO)  reco-ds.  Further  editing,  quality  control 
«»c.  re  samp  ing  reduced  the  iata  to  the  event  libraries  cataloged  in  Tables 

1,  2  and  3.  Sample  intervals  for  the  edited  data  are  144  msec  (Grapevine) 
and  72  msec  (UBO). 

Master  record  numbers  refer  to  numbers  assigned  to  the 
data  samples  edited  from  field  tapes  and  transferred  to  TIAC  tapes  in  the 
first  stage  of  editing.  Grapevine  data  are  stored  on  TIAC  tapes  1453.  1293 
and  1 14  under  various  formats.  Table  1  lists  the  28  events  from  the 
Grapevine  library  which  were  selected  for  use  in  the  analysis  program. 
Additional  library  data  not  used  in  the  analysis  are  listed  in  Table  2. 

From  the  UBO  field  tapes.  15  events  were  selected  for  analysis  and  written 
on  tape  1981. 
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Formats  of  the  library  data  are  as  follows: 


I 


1 

I 

I 

0 


n 
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•  Tape  1453 

(1)  18-trace  records  containing  low-gain  recordings 
of  the  surface  vertical,  six  deep  verticals,  and  two 
horizontals  in  traces  1  through  9  and  high-gain  re¬ 
cordings  in  traces  19  through  18  (eight  events). 

(2)  21 -trace  records  containing  outputs  of  the  surface 
vertical  and  six  deep  verticals  in  traces  1  through 
7.  Three  events  were  stored  in  each  record  simi¬ 
larly  assigning  traces  8  through  14  to  a  second 
-,-ent  and  traces  15  through  21  to  a  third  event 

(46  events). 

•  Tape  1293 

7 -trace  records  containing  whitened  outputs  of  the  sur¬ 
face  vertical  and  six  deep  vertical  in  traces  1  through  7, 
Some  14-records  containing  a  second  event  in  traces 
8  through  14. 

•  Tape  114 

Same  as  Tape  1293  except  no  record  contains  more 
than  one  event. 

•  -'ape  1981 

29-trace  records  containing  10  on-line  processor 
outputs  in  traces  1  through  10,  10  shallow-buried 
verticals  in  traces  11  through  20,  surface  vertical 
2-10  in  trace  21,  two  horizontals  in  traces  22  and  23 
and  six  deep-well  verticals  in  traces  24  through  29. 
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SECTION  II 

multichannel  filter  design  for  grapevine 


For  the  Grapevine  installation,  6-channel  Wiener  filters 
were  designed  to  operate  upon  the  outputs  of  the  aix  deep  instruments  to 
optimally  estimate  the  signal  that  would  be  observed  by  a  surface  seis¬ 
mometer  in  the  absence  of  noise.  The  models  used  to  derive  the  filter 
responses  are  listed  in  Table  4. 


Table  4 


multichannel  filters  for  grapevine 


Filter  No. 

Noise  Model 

Signal  Model 

1 

Average  of  12  noise  samples 

Honshu  earthquake, 
with  surface  trace 
replaced  by  time- 
shifted  output  of 
velocity -filtering 
the  six  deep  out¬ 
put  traces 

2 

Average  of  12  noise  samples 

Honshu  earthquake, 
with  actual  surface 
trace  used  as  de¬ 
sired  signal 

3 

Average  of  6  intense 
noise  samples:  G1 -G6 

Theoretical  signal 

4 

Average  of  6  quiet 
noise  samples:  G7 -G12 

Theoretical  signal 

5 

Average  of  12  noise 
samples 

Theoretical  signal 

13 
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For  each  record,  relative  instrument  gain  corrections  were 
applied  based  on  calibration  signals.  Then  a  whitening  filter  was  designed 
from  the  autocorrelation  of  the  surface  trace  and  applied  to  all  traces  in  the 
record.  The  whitened  data  were  used  to  compute  experimental  signal  and 
noise  correlation  sets  for  input  to  the  MCF  design  programs.  In  addition. 

a  set  of  signal  correlations  was  derived  theoretically  for  a  vertically  incident 
P-wave  signal. 

Outputs  of  three  of  the  filter  sets  are  displayed  in  Figures  3 
through  H  along  with  playbacks  of  the  original  data  (after  the  deconvolution 
process  described  above).  Since  the  three  filter  sets  designed  from  theo¬ 
retical  signal  correlations  had  almost  identical  outputs,  only  one  of  them 
is  presented  in  this  report.  For  the  Grapevine  data,  all  design  and  appli¬ 
cation  of  multichannel  filters  was  performed  on  the  TLAC  computer. 

At  Grapevine,  a  severe  attenuation  of  seismic  noise  with 
depth  of  burial  is  observed.  The  variability  of  the  attenuation  is  illustrated 
by  the  obvious  difference  between  noise  sample  552  (Figure  3)  and  573 
(Figure  4)  which  were  recorded  at  0609  and  1057,  respectively,  local 
standard  time.  In  a  general  wideband  sense,  the  improvement  in  signal- 
to-noise  ratio,  obtained  by  applying  multichannel  filtering  to  vertical  array 
outputs  was  not  significantly  greater  than  the  improvement  obtained  by 
burying  a  single  instrument  (Figures  3  through  6).  Figure  7  is  a  repro¬ 
duction  of  Figure  5  with  analog  playback  filters  applied  to  reject  energy 
below  1,66  cps.  In  the  high-frequency  range,  Filter  2  shows  a  signal- 
to-noise  improvement  of  roughly  10  db  relative  to  the  quietest  deep-well 
seismometer.  These  observations  suggest  that  vertical  arrays  offer 
promise  of  some  useful  improvement  above  1.5  cps. 
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Figure  7.  Low-Cut  Filtered  Playbacks  of  Alaska  Earthquake 
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Figure  8. 


Synthetic  Teleseism  and  the  Responses  of  Three 
Filter  Sets 


Grapevine 


Filters  1  «md  2,  designed  from  experimental  signal  statistics, 
produce  a  serious  signal  distortion.  This  effect  is  shown  by  the  results  of 
applying  the  filters  to  a  synthetic  telesism  (Figure  8).  Since  very  noisy 
records  were  used  as  estimates  of  the  desired  signal  in  designing  these 
two  filters  it  is  to  be  expected  that  such  distortion  should  appear.  Con¬ 
versely,  these  filters  provided  some  signal-to-noise  improvement’  whereas, 
signal-to-noise  ratios  ir  the  outputs  of  Filters  3,  4  and  5  were  generally 
worse  than  those  in  the  outputs  of  single  deep-well  seismometers.  It 
appears  that  the  synthetic  signal  used  to  derive  these  filters  wng  different 
from  the  signals  actually  observed,  possibly  because  of  differences  in  the 
coupling  factors  between  the  deep-well  seismometers  and  the  surrounding 
formations.  Standard  calibration  techniques  are  not  sufficient  to  resolve 
this  problem  since  they  measure  the  response  of  the  system  to  motions  of 
the  seismometer  but  not  the  recponse  to  motions  of  the  surrounding  rock. 
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SECTION  III 

SIGNAL  ENHANCEMENT  FILTERS  FOR  UBO 


Recordings  of  the  Yukon  and  Mid-Indian  Rise  earthquakes 
and  six  noise  samples  were  used  to  obtain  experimental  correlation  statitics 
lor  MCF  design.  From  each  of  the  right  records,  a  9-trace  record  was 
generated  which  contained  the  outputs  of  the  3 -component  surface  insta’lation 
and  the  six  deep  seismometers.  In  each  case,  a  whitening  filter  was  de¬ 
signed  for  the  first  trace  (surface  vertical  instrument)  and  applied  to  all 
i iine  traces  in  the  record.  (For  the  noise  samples,  autocorrelations  were 
computed  over  the  entire  record;  for  the  signals,  shorter  gates  were  used.) 

The  whitened  data  were  then  used  to  compute  experimental  signal  and  noise 
correlation  sets.  In  this  instance,  entire  records  were  used  to  compute 
both  signal  and  noise  correlations. 

The  eight  sets  of  correlation  functions  were  transcribed 
to  IBM  format  tape  so  that  they  could  be  used  by  filter  design  programs 
operating  on  the  IBM  7044.  In  addition,  theoretical  signal  correlations 
were  computed  from  the  equations  for  wave  propagation  in  a  layered 
elastic  medium,  assuming  vertical  incidence  and  including  the  effect  of 
all  reflected  waves.  In  all,  eight  sets  of  multichannel  filters  were  designed, 
and  of  these,  seven  sets  were  transcribed  back  to  the  TIAC  computer  and 
applied  to  the  suite  of  15  records  listed  in  Table  3. 

Noise  statistics  for  each  filter  design  were  obtained  by 
averagin;  the  corre'ation  sets  computed  from  the  six  noise  samples.  Three 
of  the  filtei  sets  were  designed  from  experimental  signal  statistics  and 
five  were  based  on  theoretical  signal  statistics.  A  difficulty  associated 
with  the  use  of  theoretical  signal  statistics  is  that  false  gains  may  result 
if  experimental  noise  statistics  are  derived  from  unequalized  instruments. 
Equalizing  a  vertical  array  is  a  complicated  problem  since  both  signal  and 
no;se  amplitudes  vary  as  functions  of  depth  and  frequency. 
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Calibration  results  (which  were  not  used  in  designing  any 
of  the  multichannel  filters)  suggested  that  relative  gains  of  the  seven 
vertical-component  seismometers  (from  surface  to  deepest)  were  1.0, 
1.0833,  0.9175,  1.0H87,  0.8229,  1.0403,  1.0359.  Crude  examination  of 
signal  amplitudes  suggested  that  actual  instrument  gains  were  1.0,  0.82, 
0.86-  ,  0.694,  0.93,  0.885,  0.840;  i.  e.  ,  it  was  observed  that  the  largest 
signal  amplitudes  appeared  in  the  outputs  of  the  instruments  which  should 
have  had  the  lowest  gains  according  to  the  calibration  data.  This  effect 
may  be  due  to  the  inability  of  the  calibration  method  to  take  into  account 
P  'ssiole  differences  in  coupling  factors  between  the  seismometers  and 
rock  formations. 

This  dilemma  was  attacked  in  the  last  filter  design  by 
applying  gain  factors  to  the  signal  autocorrelations  so  as  to  describe  a 
signal  as  seen  by  a  set  of  instruments  with  random  differences  in  gain. 

The  second  and  third  filter  sets  were  9-channel  filters,  designed  to  operate 
on  the  outputs  of  the  nine  instruments  so  as  to  estimate  the  signal  portion  of 
the  output  of  the  surface  vertical-component  instrument.  Since  it  was 
readily  apparent  that  no  significant  advantage  was  obtained  by  including  the 
two  horizontal  components  (the  filter  responses  for  the  two  horizontal 
seismometers  were  nearly  zero),  the  other  filter  sets  were  7-channel 
filters  operating  on  only  the  vertical -component  instruments. 

Characteristics  of  the  UBO  multichannel  filter  sets  are  given 
in  7  Jble  5.  Of  the  eight  filter  sets,  only  the  last  one  provided  an  output 
which  was  superior  to  a  single  seismometer  (Figures  9  through  12).  This 
d  ng  suggests  that  the  multichannel  filter  design  for  vertical  arrays  should 
include  an  allowance  for  uncertainties  in  instrument  coupling  factors. 
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Figure  10.  Playbacks  of  Rat  Islands  Earthquake  Observed  at  UBO 


IF 


In  general,  FUter  106  has  provided  an  output  which  is 
superior  to  that  of  a  single  deep-well  seismometer,  since  it  is  possible 
to  eliminate  interference  effects  caused  by  the  surface  reflection.  The 
spectral  estimates  plotted  in  Figure  13  imply  a  noise  reduction  by  Filter 
10^  of  about  8  to  10  db  in  the  frequency  range  above  1.  5  cps,  relative  to 
a  single  surface  instrument.  On  the  other  hand,  the  attenuation  of  seismic 
noise  with  depth  is  so  small  at  UBO  that  the  filter  output  does  not  show  a 
significant  wideband  advantage  over  a  single  seismometer  at  the  surface. 
However,  Figure  14  shows  that  a  considerable  improvement  can  be  ob¬ 
tained  in  the  frequency  range  above  1.  5  cps.  It  appears  that  useful 
high-frequency  singal  enhancement  can  be  obtainable  from  a  well  calibrated 
vertical  array. 


Figure  13.  Estimated  Response  of  Filter  Set  106  to  Six  Noise  Samples 
(Relative  to  a  Single  Surface  Seismometer) 
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SECTION  IV 
NOISE  ANALYSIS 


Correlation  sets  computed  from  12  Grapevine  nni8e  samples 
and  f  1  h°  noise  8  *mPle8  were  transcribed  to  the  IBM  70<.  i  for  further 
analysis.  For  each  site,  a  set  of  averaged  correlation  funct  ons  was  used 
to  compute  squared  coherence  functions  between  all  pairs  of  vertical  seis¬ 
mometers  in  the  vertical  array.  These  coherence  functions  are  displayed 
in  Figures  15  and  16.  It  is  observed  that  instruments  separated  by  1000  ft 
display  coherence  greater  than  0.  5  throughout  the  frequency  band  0  to  2  cps. 
The  coherence  decreases  with  increasing  frequency  above  2  cps.  As  the 
separation  is  increased,  the  frequency  band  containing  coherent  noise  energy 
becomes  smaller.  No  useful  coherence  is  found  between  the  surface  instru¬ 
ment  and  the  deepest  instrument. 

Attempts  were  made  to  perform  a  modal  analysis  of  the  noise 
fie.is  by  fitting  combinations  of  theoretical  multichannel  correlation  sets 
to  the  experimentally  derived  multichannel  noise  correlation  sets.  As  has 
been  described  in  Semiannual  Report  No.  4.  theoretical  sets  were  derived 
for  models  with  white  spectra  as  observed  by  a  surface  seismometer.  The 
"white"  correlation  sets  were  filtered  (one  filter/correlation  set)  with 
appropriate  filters  so  that  they  could  be  combined  to  produce  the  best  match 
to  the  experimental  correlation  sets.  Filter  amplitude  responses  were 
interpreted  as  estimates  of  the  power  spectra  of  the  corresponding  modes 
contributions  to  the  averaged  whitened  experimental  noise  models. 

Figures  17  and  18  show  theoretical  correlation  sets  derived 
for  Grapevine  and  UBO.  In  each  case,  theoretical  models  included  verti¬ 
cally  incident  P-waves;  surface-wave  modes  0,  1  and  2;  and  incoherent 
noise  with  equal  power  at  all  depths.  In  addition,  in  the  Grapevine  analysis 

only,  a  sixth  model  of  noise  was  incoherent  energy  appearing  only  at  the 
surface. 
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Ft  r  h  “**’  h‘ee  se“'  ot  filters  were  designed  to  operate 
he  theoretical  correlation  sets.  Figure  1,  shows  the  of  eJ[perimental 

eo.se  correlations  from  Grapevine,  the  optimum  estimates  obtained  by 
multichannel  filtering  the  six  theoretical  correlation  sets  with  21-  31- 

-d  4, -p0i„t  mters  and  ,he  est.raation  ._rrora  Figur<=  2o  Bh(>w>  (he 

of  experimental  noise  correlations  from  UBO.  the  optimum  estimates 
Obtained  by  multichannel  filtering  the  five  theoretical  correlation  set,  with 
31.  and  49 -point  filters,  and  the  corresponding  errors.  The  noise 
correlations  in  Figure  20  were  produced  by  applying  gains  to  the  origins, 
correlation  sets  to  compensate  for  instrument  responses  estimated  visually. 

The  filter  responses  have  been  used  as  estimates  of  the 
spectra  contribution,  of  the  various  modes.  The  results  for  Grapevine 
are  shown  in  Figure  21.  All  three  filter  sets  indicate  that  the  seismic 
noise  a,  the  surface  is  dominated  by  incoherent  energy  which  doe.  not  appear 
epth.  The  most  important  coherent  mode,  are  found  to  be  P-wave.  below 
•  cps.  mode  i  surface  wave,  between  0.  5  and  1  cp.  and  mode  0  surface 
wave,  above  1  cps.  Since  Grapevine  is  a  noisy  site  which  suffers  from 
ocaily  generated  cultural  disturbances,  i,  doe.  no,  seem  reasonable  tha, 
he  seismic  noise  should  be  dominated  by  high-velocity  P-wave.  a,  any 

requency.  Therefore,  the  result,  obtained  suggest  that  the  method  of 
mode  separation  may  not  have  any  validity. 

Result,  obtained  for  UBO  are  shown  i„  Figure  22.  Again. 

It  is  ound  tha,  P-waves  appear  to  be  the  dominant  noise  propagation  mode 
at  low  frequencies. 

The  process  of  mode  separation  depenos  upon  recognizing 
epth-dependence  functions  corresponding  to  the  theoretical  models.  Since 
differences  between  depth-dependence  functions  for  two  different  modes  can 

/  suotle,  it  ,s  to  be  expected  that  calibration  error,  could  produce  serious 
errors  xn  a  mode  separation  analysis. 
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To  test  this  hypothesis,  the  UBO  mode  sepa.-tion  procedure 
was  repeated  with  the  experimental  correlations  gain-corrected  according 
to  the  values  derived  from  visual  measurement  of  teleseismic  signal  ampli¬ 
tudes.  Since  there  was  a  great  difference  between  the  two  seti  of  estimates 
of  instrument  responses,  it  is  to  be  expected  that  the  mode  separation 
results  should  be  altered  drastically.  However,  the  results  plotted  in 
Figure  23  show  that  the  method  is  virtually  insensitive  io  large  errors  in 
instrument  gain.  This  finding  also  casts  doubt  m  the  method  of  analysis. 


figure  15.  Iwo-Channel  Coherence  Functions  Obtained  by  Averaging 
Noise  Samples  from  Grapevine 
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Figure  21.  Estimated  Contributions  of  Propagation  Modes  to  Grapevine  Noise: 
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Figure  22.  Estimated  Contributions  of  Propagation  Modes  to  UBO  Noise 
Corrected  for  Visually  Estimated  Instrument  Responses:  (a) 
21-Point  (1.44  sec)  Filters,  (b)  31-Point  (2.  16  sec)  Filters, 
(c)  49-Point  (3.456  sec)  Filters 
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Figure  23, 


Corrected  f  *  “  °f  ProPagati°n  Modes  to  UBO  Noise 

Instrument  Gains  Derived  from  Calibration  Sig- 

mte  !  *44  Sec)  Filters'  (b)  31-Point  (2.  16  sec) 

Filters,  (c)  49-Point  (3.456  sec)  Filters 
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